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Abstract. We investigate the effects of aperiodic interactions on the 
critical behavior of an interacting two-polymer model on hierarchical lattices 
(equivalent to the Migadal-Kadanoff approximation for the model on Bravais 
lattices), via renormalization-group and tranfer-matrix calculations. The exact 
renormalization-group recursion relations always present a symmetric fixed point, 
associated with the critical behavior of the underlying uniform model. If the 
aperiodic interactions, defined by s ubstitution rules, lead to relevant geometric 
fluctuations, this fixed point becomes fully unstable, giving rise to novel attractors 
of different nature. We present an explicit example in which this new attractor 
is a two-cycle, with critical indices different from the uniform model. In case of 
the four-letter Rudin-Shapiro substitution rule, we find a surprising closed curve 
whose points are attractors of period two, associated with a marginal operator. 
Nevertheless, a scaling analysis indicates that this attractor may lead to a new 
critical universality class. In order to provide an independent confirmation of the 
scaling results, we turn to a direct thermodynamic calculation of the specific-heat 
exponent. The thermodynamic free energy is obtained from a transfer matrix 
formalism, which had been previously introduced for spin systems, and is now 
extended to the two-polymer model with aperiodic interactions. 
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1. Introduction 

In some recent publications ^-0], we have used renormalization-group (RG) and 
transfer-matrix (TM) techniques to investigate the effects of aperiodically distributed 
(but not disordered) interactions on the critical behavior of ferromagnetic spin models. 
In a real-space renormalization calculation for simple hierarchical structures, we have 
written exact recursion relations in order to show that relevant geometric (aperiodic) 
fluctuations play a very similar role to disorder. These calculations lead to the 
formulation of an exact extension for deterministic, aperiodic interactions of the well- 
known Harris criterion for the relevance of disorder in a form coincident with 
Luck's heuristic, general derivation of this extension Also, we have shown that 
relevant geometric fluctuations give rise to distinct critical exponents, associated with 
the appearance of new attractors in parameter space. The independent transfer- 
matrix calculations have confirmed these results and provided deeper insight on more 
refined details of the thermodynamics of aperiodic spin systems, such as log-periodic 
oscillations of thermodynamic functions. 

Now we consider a model of two directed polymers on a diamond hierarchical 
lattice, with an aperiodic, layered distribution of interaction energies, according to 
various substitution rules (see IH| for extensive reviews of substitution sequences 
and their applications to statistical models). Although the qualitative description of 
the critical behavior is essentially similar, the exact renormalization-group recursion 
relations in parameter space turn out to be much simpler as compared to the 
calculations for spin systems. The case of two-letter substitution rules is particularly 
simple. For irrelevant geometric fluctuations, the critical behavior is governed by a 
symmetric fixed point, with no changes with respect to the uniform case. For relevant 
fluctuations, we show that this symmetric fixed point becomes fully unstable, and 
the critical behavior is associated with a novel two-cycle attractor of saddle-point 
character. However, for more complex substitution rules, such as the Rudin-Shapiro 
sequence of four letters (which is known to mimic Gaussian random fluctuations in a 
sense [H]) there appears a surprisingly rich structure in the four-dimensional parameter 
space. Besides the expected symmetric fixed point, there are non-diagonal fixed points 
and continuous lines of two-cycle attractors, associated with a marginal operator, 
which might give rise to non-universal critical exponents. A scaling analysis indicates 
that this structure is responsible for a novel critical universality class. In order to 
check these results, we resorted to an independent thermodynamic calculation, on the 
basis of a transfer matrix scheme. 

The layout of this paper is as follows. In Section 2, we define the polymer model, 
and present some renormalization-group calculations for two- and four-letter (Rudin- 
Shapiro) substitution rules. We show the existence of some surprising structures in 
parameter space, and, in particular, discuss a number of scaling results for the critical 
behavior. We then proceed to formulate an extension of the transfer matrix scheme for 
a two-polymer model. Although this technique has already been used for spin systems, 
its extension for this new situation requires a considerable amount of analytical work, 
as described in Sections 4 and 5. In Section 3, we just formulate the transfer-matrix 
scheme for a two-polymer model. In Section 4, we study in great detail the algebraic 
structure of the transfer matrices, in order to unveil, in Section 5, the existence of 
a recursion relation for the eigenvalues of the transfer matrices associated with two 
successive generations of the hierarchical structure. In Section 6, we use the recursion 
relations in order to write down explicit thermodynamic functions. In Section 7 we 
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present the results for aperiodic models, which are compared to the renormalization- 
group predictions. 

2. The interacting polymer model 

The binding-unbinding phase transition in a disordered model of two directed and 
interacting polymers on a hierarchical lattice has been investigated by Mukherji and 
Bhattacharjee [S], and we follow these authors on the definition of the model. We 
simply place two directed polymers on a diamond hierarchical. They start at one end 
of the lattice and stretch continuously to the other end (the root sites). There is an 
attractive interaction, — e, whenever a bond of the lattice is shared by a monomer 
of each polymer. This energy can be made to depend on the position of the bond 
along a branch, in a random or deterministic fashion. Note that, although seemingly 
artificial, the model is nothing else than a bona-fide Migdal-Kadanoff approximation 
for the same interacting problem on a genuine Bravais lattice. 

In the basic cell of a diamond lattice, with q branches and p bonds per branch, 
there are configurations of energy — pe, where the two polymers occupy the same p 
bonds of a branch, and configurations of zero energy, where the polymers stretch along 
different branches. Using the Boltzmann factor y = exp (/3e) and the combinatorial 
coefficient — q{q — 1) /2, it is easy to write the RG recursion relation 

y' = V + ^. (1) 
q q 

Taking p — 2, for example, we see that besides the trivial fixed points, y* — 1 and 
oo, associated with zero and infinite temperatures, there is a nontrivial fixed point, 
y* = q — 1, which is physically acceptable for q > 2 (there is no phase transition on 
the simple diamond lattice with q = 2 branches). Also, from the linearization of the 
recursion relation about the nontrivial fixed point, we have the thermal eigenvalue 

A=^. (2) 

In order to obtain the specific-heat exponent, a, one should note that, as the 
polymers are one-dimensional objects, the thermodynamic extensivity of this model 
relates to the polymer length, instead of the volume. Thus, the important density is the 
free energy per monomer, which is assumed to behave according to the fundamental 
scaling relation 

fiy')=pf{y)- (3) 

From this relation, we have the critical exponent associated with the specific heat. 

This quantity will be used to characterize the possible universality classes of the model, 
depending on q, p and the presence of aperiodically distributed interactions. 

Consider, for example, for p = 2, a layered distribution of interactions 1 , and 
Cf,, chosen according to the two- letter period-doubling sequence, 

a — > ab, b — > aa. (5) 

In Figure 1 we give an example of the construction of a simple diamond lattice with 
this kind of aperiodicity, starting from the letter a. Along each branch the interaction 
energies are distributed according to the letters of the aperiodic sequence generated 
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Figure 1. Initial stages of the construction of a Migdal-Kadanoff hierarchical 
lattice with q = 2 branches and p = 2 bonds per branch, and layered aperiodic 
interactions according to the period doubling rule, a — > ab, b ^ aa (letters a and 
b indicate the two possible values of the interaction energy, and ej). 



by the recursive application of the rule. The same arguments as used in the last 
paragraph to derive Eq. for the uniform case, lead to a pair of recursion relations, 

Va^-VaVb^ , (6) 



q 



and 



where = exp(/3ea,6), and ea,b > is the interaction energy at bonds of types a 
and b, respectively. For Ua = Vb — y, we recover the recursion relation associated with 
the uniform model. 

It is easy to see that, for q > 2, there is no physically acceptable nontrivial fixed 
point, except the symmetric fixed point, yab ~ V* = 9 ~ 1- The linearization of the 
recursion relations © and (Q in the neighborhood of this symmetric fixed point leads 
to the matrix form 

with eigenvalues Ai — 2y* /q — 2 {q ~ 1) /q, and A2 = —y*/q = —{q— l)/?- Therefore, 
if q > 2, we have Ai > 1 and IA2I < 1, which shows that the geometric fluctuations are 
completely irrelevant in this case (since Ai is the same eigenvalue associated with the 
non-trivial fixed point of the underlying uniform model, and the behavior is irrelevant 
along the other direction). 

We now turn to a more interesting case. Consider a diamond lattice with p = 3 
bonds and q branches. Suppose that the (layered) interactions are chosen according 
to a period-3 two-letter sequence, a — + abb and b — > aaa. The new recursion relations 
are given by 

y'a = -VaV! + (9) 



and 



y'b = -y'a + ^- (10) 



q q 

For q > 3, there is just a single nontrivial fixed point, at a symmetric location, 

1 1 



yl^yt^y* ^-^ + ^V^q^- (u) 
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The Unearization in the neighborhood of this fixed point leads to the matrix equation 

\2 



Ayb 



with eigenvalues 



q 2q 



2q-l- ^iq-3 



and 



*2 



A2 = -2 



y 
q 



1 r 

q 



2(7 - 1 - v/49-3 



(12) 



(13) 



(14) 



For 3 < q < 3 + a/5, it is easy to see that Ai > 1, and IA2I < 1. As in the 
case of the simple diamond lattice with p — 2 bonds, geometric fluctuations are 
irrelevant and the critical behavior is identical to the uniform case. However, for 
g > 3 + a/S = 5.236068..., we have IA2I > 1, and the symmetric fixed point becomes 
fully unstable, and, therefore, cannot be reached from arbitrary initial conditions (the 
values of ea and et). For example, for g = 5, we have y* = = y* = 1.561552..., 
with eigenvalues Ai = 2.140568... and A2 = 0.951363... < 1. For q — 6, however, 
we have y* = y* ^ y* = 1.791287..., with eigenvalues Ai = 2.573958... and 
A2 — 1.143981... > 1. But, as in the case of spin models on hierarchical lattices 
, there is a two-cycle in parameter space. It is easy to numerically locate this cycle 
at {yl,y;) = (1.419001..., 2.267305...) and (2.049103..., 1.309541...), with eigenvalues of 
the linearized second iterate given by Ai = 2.624300... > 1 and A2 = 0.772598... < 1. 
A new critical universality class is therefore expected to be defined by this attractor. 

The behavior in parameter space is much more interesting if we consider the 
Rudin-Shapiro, four-letter substitution rule, a ac, b — > dc, c ~^ ah, d — > db. 
Consider a simple diamond lattice with p — 2 bonds and q branches. It is 
straightforward to write four algebraic recursion relations. 



y'a 



yayc/q + (g - 1) /g. 



y'b 



ydyc/q + (q- l) /q, 



y'c = yayb/q + (q - l) /q, y'd = ydyb/q + (g - 1) /g, 

which lead to the symmetric fixed point 

y*a = yt = yc ^ y*d y* q - ^■ 

From the linearization about this fixed point, we have the eigenvalues 



Ai 



1 



A, 



V2 



q-i 



A.. 



-x/2 



q-l 



A4 0. 



(15) 
(16) 



(17) 



(18) 



q q q 

The introduction of aperiodic interactions becomes relevant for the simple diamond 
lattice if g > 2 -I- \/2 = 3.41..., which corresponds to IA2I = IA3I > 1. 

Recursion relations 1)15(1 and ((16|l are so simple that we can perform a number 
of detailed calculations. In particular, it is easy to show the existence of additional, 
non-diagonal fixed points, given by 

1 



4(g-l) 
q-1 

q — — 

Vd^a 



q (g2 - 2g + 2) ± q^ [q' ^ 2) (q^-Aq + 2) 



(19) 
(20) 
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The Jacobian matrix associated with the Unear form in the neighborhood of these 
fixed points can be written as 

/ y* \ 

vt • ^^'> 
V y*, yi I 

Besides two trivial eigenvalues, A3 = and A4 = —1, there is an additional pair of 
eigenvalues (|Ai| > 1 and IA2I < 1) given by the solutions of the quadratic equation 

2q^K^ - V {q^ + q-l) I^ + q^ + 8q^ - 4q^ - 8q + 4 = 0. (22) 

As IA4I = 1, we have a typical case of marginal behavior, which cannot be analysed 
without resorting to higher-order calculations. However, the marginal operator does 
not give rise to a line of fixed points, as could be expected. Instead, with an additional 
algebraic effort, it is possible to show the existence of a continuous line whose points 
are two-cycles, by solving the polynomial equations 

ya,b,c,d = ya,b,c,d (ya, yb, yc, yd) , (23) 

where y" is the second iterate of the recursion relations. Given any q > 2 + a/2, these 
equations lead to a pair of one-parameter algebraic curves, which meet smoothly and 
form a single, non-intersecting closed curve, containing the non-diagonal fixed points. 
Any point belonging to this closed curve is mapped into another point on the curve 
upon one iteration of the recursion relations, and back to itself upon a further one. We 
have used algebraic computation to check this result very thoroughly. As an example, 
for q = 4, we have the equations 

.. = 3^, (24) 

_ 70t^ + 71t-lU±g{t) 
•^t[Qt^-im + 7S>±g{t)y 

6t^-128/ + 114±g(t) 

wta) — 

yd = t, (27) 

g{t) = {36t* - 222Qf^ + 15529*2 - 27132t + 12296)^/^ (28) 

and the parameter t taking values between 1.5142 . . . and 5.4112 In Figure 2 we 

show a three-dimensional projection of this attractor. The linearization of the second 
iterates about any point of the curve leads to the eigenvalues 0,-1, and a conjugate 
pair, I All > 1 and IA2I < 1 (the values of which do not depend on the point about 
which linearization is being carried out). The effects of the marginal eigenvalues on the 
specific- heat exponent, related to the existence of an extended attractor in parameter 
space, have to be checked very carefully, so we turn to the direct thermodynamic 
analysis of the free-energy singularity. 

Before proceeding to the transfer-matrix calculations, however, it is worth 
remarking that, from a broad renormalization-group perspective, relevance and 
irrelevance of aperiodic distributions of couplings is related to the existence or not 



and 



with 
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Figure 2. Three-dimensional projection of the attractor of the RG recursion 
relations for the Rudin-Shapiro aperiodic sequence, in a lattice with p = 2 and 
g = 4. The non-diagonal fixed points are shown as circles. 



of a second eigenvalue with modulus larger than unity. The general structure of the 
recursion relations can be used to derive a relevance criterion which is ultimately based 
on the geometry of the lattice (that is, the values of p and q), and some measure of 
the "strength" of the aperiodic fluctuations, such as the wandering exponent [H]- 

3. Transfer matrix formulation 

One of us has successfully used a transfer matrix (TM) technique to obtain the 
thermodynamic properties of several spin models on fractal lattices 3 , 4 . The essential 
step of this scheme consists in the derivation of maps relating the eigenvalues of the 
transfer matrices associated with two subsequent generations, G and G + 1. In a 
certain sense, this formalism is equivalent to a method used by Derrida et al. |10| to 
establish a map for the free energy, although it enables the calculation of a correlation 
length, which turns out to be very useful to locate the critical temperature (in spin 
systems). 

In order to introduce the transfer matrix formulation, let us define the model 
in a mathematically more precise way. Again, we will consider a simple diamond 
hierarchical lattice (remembering that it has q branches and p bonds per branch in 
each basic cell). At a generation G of the hierarchical construction, the lattice has 
q'^ branches, each of them formed by p'^ bonds. A polymer A, extending from one 
root site to the other, is formed by p^ connected monomers, and represented by a 
directed continuous path between the two root sites. A given monomer, labeled i 
(0 < z < p*^ — 1), occupies one of the q^ available branches at the zth position along 
the path. We may define a numbering for the q'~^ branches of the lattice, and let the 
variables Oi indicate which one is occupied by the ith monomer; clearly, 1 < Oi < q^ , 
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Figure 3. Periodic chain of liicrarcliical cells joined by the root sites, each 
consisting on the first generation of a diamond lattice with p = 2 and q = 2. This 
is the initial stage of the construction of the transfer-matrix formalism for the 
interacting polymer model. 



with analogous definitions for the other polymer, B. The two polymers interact at 
position i, with energy — < 0, if the ith monomers of the two distinct polymers 
occupy the same bond (note that the energy depends only on the position i, since we 
will consider layered interactions). If the zth monomers occupy distinct bonds, the 
interaction energy vanishes. 

This definition of the two-polymer model can be summarized by the Hamiltonian 



where S (a^, bi) indicates a Kronecker delta. 

In Figure 3, we illustrate the simplest case, for G = 1, p = 2 and q = 2. Note 
that, in fact, we are considering a periodic chain of N hierarchical cells, each one 
grown up to generation G. For a single cell, as each monomer can independently 
occupy any of the q branches, there are possible configurations for each position i 
(remember that, as p = 2, i can have only two values in each cell), so that a total of 
16 possible states can be devised for this specific situation. However, we have to take 
into account that each polymer is required to form a continuous path between the root 
sites, so that several configurations have to be excluded. Formally, we may calculate 
the partition function including these forbidden configurations, if we introduce in the 
Hamiltonian an additional term of the form 77 Vi(ai, a j+i, 6j, 6j+i), with rj ^ 00, 
such that Vi = for the acceptable configurations, and > whenever Oj, aj+i and 
bi, bi+i are not properly constrained. The explicit form of this potential is somewhat 
cumbersome, and will not be given here. 

Although it is straightforward to write a partition function for the particular case 
illustrated in Figure 3, 



the calculation of Zq, for arbitrary values of G, and hence of the thermodynamic 
properties of the model, represents a much more difficult task. This is the reason to 
invoke the transfer matrix technique. However, this problem of interacting polymers, 
with interaction energies depending on monomer positions along the bonds, requires 
a completely new definition of the transfer matrices. 

For the sake of simplicity, we restrict the presentation of the formalism to the 
homogeneous (uniform) system, which requires a smaller number of different types of 
elementary TM's. It is straightforward to work out an extension for the more complex 
situation of a model with aperiodic interactions. 




(29) 



Zi = {2e'^''°e^'' + 2) 



N 



(30) 
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With the inclusion of the infinite-energy term, Eq. H29|l can be written in the 
symmetrized form, 



p 

e 



T-Lg = --^ X [(5(a^,6i) + (5(aj+i,6i+i)] +77 ^ V^{a^,ai+i,bi,bi+i) , (31) 

1=0 j=0 

where we impose periodic boundary conditions, = oq and b^c = 60. Now we 
define the g^*^ x q^'^ matrices, 

(M^'')aib,;a,+ib,+i = °^p|y ['^("»'^'') + 5{a.i+i,bi+{)\ +/3?7T^i(aj,ai+i,6i,6i+i)| , (32) 

and recall that the index i ranges from to p*^ — 1 (the length of each polymer in 
generation G is p*^ ). Also, it should be remarked that, in the double indices of M^^-*, 
each term Ojbi must be interpreted as a tensor product of the variables and 6^, each 
of them taking q'^ independent values. The Boltzmann weights at each M^^' relate to 
the attracting energy between two monomers and to whether the ith monomer, placed 
at the bond Oi (bi), can be linked to the (i + l)th monomer at the bond ai+i(6i+i), 
without violating the continuity constraint. The definition of the TM's in Eq. (|32|l 
leads to the formal identification of the partition function with the trace of a product 
of TM's, 

Zg = Tt Y[ ^''g = TrMc. (33) 

i=Q 

Note that in this equation we considered = 1, since periodic boundary conditions 
have already been enforced. 

In the following we will show that M^j has just one non-zero eigenvalue, which 
we call riG, and which is obviously its trace. For the sake of clarity, we state now the 
result that rjG is given by a recursion relation, 

VG = q [vh-i + iq- l)x%^i] , (34) 

with 

XG = g(^'-^), (35) 
and rjo — y^. The next two sections are devoted to the derivation of this result, and 
some readers may be interested on skipping directly to Section 6, where we establish 
recursion relations for the relevant thermodynamic functions. 

4. Detailed structure of the matrices 

The essential difficulty in the evaluation of Zg has been shifted into the calculation of 
the eigenvalues of a product of different kinds of TM's, . In order to emphasize the 
main ideas of the method, let us explore in greater detail the structure of the matrices 
in the case illustrated in Figure 3, where now the partition function is just the one 
given by Eq. (|30|l . with ei — 62 (remember we are considering uniform interactions for 



the time being). Defining now y — exp (/3e/2) , we note that the TM's M''*'' 
the two distinct forms. 



G=l 





y' 


y 




y\ 




y 


1 


y 


1 




y' 


y 


y" 


y 


I 


y 


1 


y 


1 / 



j M ^ y y\ =t (i)^ t(i) 



I ^L^^Jr, (36) 
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Figure 4. Second stage of the construction of the transfer-matrix formaUsm. 
Now, each cell on the periodic chain is a diamond lattice of the second generation. 



and 



/ \ 

10 



\ 1 / 



1 
1 



(37) 



depending on whether the bonds at sites i and i + 1 meet at single vertex, where all 
4 bonds are connected, or at two vertices, where the bonds are connected pairwise. 
Comparing Eqs. H36() and H37|) . we note that several matrix elements in Eq. H36|l . 
which are equal to y and 1, are replaced by O's in Eq. They result from the 

presence of the term exp(— /Jry) ^ in the Boltzmann weights, indicating that the 
polymer that arrives at a vertex where only two bonds meet cannot jump to a bond 
not incident to that vertex. For this simple situation, Eq. H33|l reduces to 



Zg=i = Tr (M^^^M^°^ 



Tr L^'L 



(i)t(o) 



(38) 



The four eigenvalues of Mi = M^^^M^^"* are given by the independent products of the 



eigenvalues of 'LY'li'^'' (which are 2 and 0) and those of JY'jf^ (which are y'^ + 1 and 
0). Thus, mJ^'mJ°^ 



^^^T.(o) (-ftrhich are 2 and 0) and those of Ji^^'^i v,"^^^^" ""^^ 

has one single non- vanishing eigenvalue, which is A = 2 (j/^ + l) . 
Restricting the analysis to p = 2, we can show that a similar result is valid for any 
g, as M^^' and mJ°^ are expressed in terms of Kronecker products of matrices with 
the same structure as those which are present in Eqs. H36|l and (|37() . All elements of 
the q X q matrix L^^^ are equal to unity; (Ji^'')i,i — y"^, (J^^^)i,; — (Ji^^)/a — y for 
I — 2,. . . ,g, while all other elements are set to unity; L^''^ is the identity q x q matrix, 

while is a diagonal matrix with y"^ in the first entry, and 1 along the rest of the 
diagonal. The only non- vanishing eigenvalue is A = g (y^ + q — 1) , which confirms the 
previous result for g = 2. 

If we go into the next generation, G = 2 (see Figure 4), the basic cell of length 2^ 

is associated with three different types of matrices, M2"'-', j = 0,1,2. They describe, 
respectively, the situations where the bonds at position i meet with those at position 
i + 1 at 1, 2 or 4 distinct vertices, each one with connectivity 2(7^, 2q and 2q^ . For 
q = 2, these matrices are 



M 



(2) 











! 




' 1 


n 



y 
y 

\ y 



y\ 
1 
1 

1 / 



(39) 
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^^^2 - 1 1 



1 1 
1 1 



and 



M 



(0) 

















n 



/ 




y 








\ 




y 


1 


















1 


1 




V 








1 


1 


/ 


/ 


y' 











\ 







1 


















1 







V 











1 





(40) 



(41) 



The presence of O's in Eqs. H4UI) and (|41|l . and the form of the matrices for larger values 
of q follow from the same kind of arguments already used to discuss the generation 
G = 1. To evaluate the partition function Z2 it is necessary to identify the order in 
which the factors M^^\ M^^^ and M^^ are multiplied to form the matrix M2. This 
can be easily realized if we recall the association of the different matrix types with the 
various kinds of vertices along the hierarchical lattice. We thus have 



M2 = m'"'m(^)m(°)m(^\ 



(42) 



from which it is straightforward to calculate Z2. 

Let us now obtain the structure of the matrices Mg for any G. This follows 
from Eq. H33|l . from the hierarchical structure of the lattice, and from the detailed 



discussion of the form of matrices M^*-* and M2''' . The first important property related 
to the structure of the lattice is that, along each branch, there are sites with different 
connectivities, and they appear according to a well-defined sequence. The difference in 
connectivities results in local variations of degrees of freedom, since the polymers can 
choose among different numbers of branches. The inner sites can be of G types, the 
connectivities of which are 2, 2q, 2q^, . . . , 2q'^~^, while the root sites have connectivity 
2q'^ . We identify the type of a site by the variable s, such that the connectivity of 
a particular site is 2q'^. Let Sg be the sequence of p*^ numbers that identifies the 
order in which the several kinds of sites appear along a branch in generation G. For 
instance, 5*1 = {0, 1} and 5*2 = {0, 1, 0, 2}. We note that Sg can be obtained by the 
concatenation of two sequences Sg-i, replacing the last symbol (G — 1) by G. Also, 
we observe that Sg contains 2*^^^ symbols 0, 2*^^^ symbols 1, and so forth, until 
one single symbol G — 1 (at the central position) and one symbol G at the rightmost 
position. As Eq. H42|) suggests, the matrix associated with a site of type s is M^^', so 
that M(3 decomposes into a product of elementary matrices in a well defined order, 



(0 



(43) 

We now study the structure of the q^'^ x q^'^ matrices Mg-* . Each of them can be 
expressed by the Kronecker product of two q'^ x g*^ matrices, and Jg, so that 



Mr 



(0)t_(1)t.(0)t.(2)t.(0)t.(1)t.(0)t(3) J^Wj^iG) 



T \^JT i>J-;t t,u;r \-^)r ^^JT t,-i-;r \^)r 
^G ^G ^G ^G ^G ^G ^ 

■t(0)t(1)t(0)t(2),(0),(1),(0) 



G 

(3) 



G "^G 



G 
(G) 



(44) 



'G "G 

= Qg «) Rg, 

where the sequence Sg controls the numbering order in both Qg and Rg. Then, we 
notice from Eqs. H39I41|I that each L^-* can be further expressed by the Kronecker 
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product of G matrices of order q, each of which is either the unit q x q matrix I or the 
constant q x q matrix K, with all elements set to unity. 

The factor K expresses allowed transitions of the polymer from a given branch 
to neighboring ones, while I indicates a restriction for the polymer to change from 

(G) 

a branch to others. So, it is easy to see that , which describes the site with 

connectivity 2q'^, is formed by the product of G matrices K, while Lq-*, related to 
sites with connectivity 2g°, is formed by products of matrices I only. If we define A*^*^ 
to be the Kronecker product of G matrices A, then we can write the general form of 



the matrices Lg-* as 

t(9) 



I®(G-g) ^ j^e 



(45) 



Also, it should be noted that, for g < G, the matrices L^-* and Lj^^^ are related by 



(9) 



' (9) 



(i) 

The matrices J W cannot be further decomposed in terms of Kronecker products. 



1(g) L 



(9) 

G-l- 



(46) 



but they can be expressed as 



t(9) _ T (9) 



H 



(9) 

G ' 



where H^' is given by 



1 



y - 



If 5 < G, it is also possible to see that 



(Hg'). 



H 



(9) 



if i = j = 1 

if^ = l, je{2,...,qS} 
if {2,...,9n , j = i 

otherwise. 
.®0, 



(47) 



(48) 



(49) 



where indicates the null q'^ ^ x q'^ ^ matrix, which appears q—1 times in the direct 
sum. This expression can also be written as 



H 



(9) 



/ 1 





V 



\ 







/ 



H 



(9) 

G-l' 



(50) 



where the first matrix is of order q. Combining this last equation with Eq. (|46|l . we 
obtain 



r(9) 



I«)L 



(9) 

G-l 



H 



(9) 

G 






V 





■(9) 

^G-1 








(9) 

-■G-l 



\ 



rig) I 

^G^l I 



(51) 





which is valid for g < G. In this expression, the O's represent null matrices of order 
9 



'-^ ^ so that Jq' is a block-diagonal matrix, with one J^lj^ block and g — 1 blocks 



LqIi in the diagonal. 
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The particular structure of the matrices and leads to a recurrence 

relation for the only non-zero eigenvalue of Mg in terms of the corresponding 
eigenvalue of Mg_i. This will be shown in the next Section. 



5. Eigenvalues of the matrix 



The eigenvalues of Mg are given by all distinct products of the eigenvalues of Qg 
and Rg- Let us first consider the eigenvalues of Qc- According to Eq. (g^, Qc 
is expressed by usual matrix products of matrices which are themselves Kronecker 
products of only two types of matrices, I and K. Then it is easy to show that 

9G-1 



(23 



1) 



G- 1, 



Qg = K (g) (g) 
Using the relation 
K" 

and the identity 

G-l 

E 

9=1 

it is possible to write Eq. (|52|l as 

/G-l 

qg = n q''-' 

The rank of K is unity, and its only one non-zero eigenvalue is q. It follows that is 
the only non-zero eigenvalue of K®*^, so that xgi the eigenvalue of Qg, is given by 



K 



®G 



(52) 
(53) 

(54) 
(55) 



XG^q^ >■ 



(56) 



Now, let us calculate the eigenvalues of Rg, defined through Eq. 14411 . First, we 
write Rg in the form 



R 



G — 



{Sg) 
G 



G 



AG) 

'g ' 



(57) 



where {Sg)i represents the /th number in the sequence Sq- Then, we note that the 



matrix J^?^ can be written as 



r(G) 



/ "^G-1 

Fg-1 
Fg-1 



Eg-1 

t(G-1) 
■^G-1 
T (G-l) 
■^G-1 



Eg-1 

t(G-1) 
^G-l 
t(G-1) 
^G-1 



V F 



G-l 



, (G-l) t(G-1) 



-•G-l 



-•G-l 



Eg-1 \ 

T (G-l) 
■^G-1 
T (G-l) 
■^G-1 



■ (G-l) 
^G-1 



(58) 



where all elements of the q'^ ^ x q^ ^ matrix Eg-i are equal to unity, with exception 
of those of the first row, which are equal to y. The matrix Fg-i is the transpose of 
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E(3_i. For g < G, we recall that Eq. uncovers the block-diagonal structure of 



3q\ so that 



2== 

n Jc 

1=1 



nlr^L^f^^' 



V 



pr2=^-l ^ (Sg), 



(59) 





with q blocks of order cp~^ in the diagonal. Let us now introduce the notation 



(60) 



with the analogous definition for the product of the matrices ■ Making use of Eqs. 
and (|5i^ . we obtain 



G-1 



G-l 



HlFg-i HLLgi-/) HlL^?-/^ 
HlFg-i HLLgr/^ HLLgL-/) 



-•G-l 



HjEg-i \ 



Let us consider the entries of this matrix. Take, for instance, 



^ t(G-1) _ 

iil1^g-i - 



n 4^-1' 



. (G-l) 
-•G-l 



l=\ 



(61) 



n 4^-1' 



n 



^G-1 



.(G-l) 



1=1 



.1=20-1+1 



(62) 



The first factor is clearly Qg-i, as the sequence Sg is equivalent to S'g-i until the 
position I = 2^~^. However, Sg is also identical to S'g-i between the positions 
I — 2*^^^ + 1 and Z = 2*^ — 1, according to the rule to generate Sq from Sg-i- Thus, 

fG— 1) 

this factor, multiplied by T-^g-i i ^-Iso leads to Qg-i, 

Qg-1- 

Following the same arguments, it is possible to show that 



TT r(G-l) 



n i^G-i) 
HjJg-i 



"-G-l' 



(63) 
(64) 



so that 



/ R2 



R 



■G — 



G-l 

HlFg-i 
HlFg-i 



HjE, 



G-l 



4G-1 

Qg 1 



HjE 



G-l 



iG-1 

Qg-1 



HjE 

Ql 



G-l 



iG-1 

Qg-1 



(65) 



G-l 



\ HlFg-i ^g-1 *°«g-i 

Now we recall that Rg is a product of the matrices J^?-* , including 3^^^ . The rank 
of this last matrix is unity, as one sees from Eq. ifsS)) . Using the Frobenius inequality 
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for the rank of matrices m, we then see that the rank of Rq is also unity. Then, the 
only non-zero eigenvalue Xq equals the trace of , so that 

XG = Trnl_, + {q-l)TvQl_,. (66) 

However, it is clear that Ti'Rq_^ = Xq_i. Also, from Eq. (|56|l . we have TrQ^_j^ = 
Xg-1, so that 

XG^Xl_, + {q-l)xl-i. (67) 

As Ao — , this equation gives rise to a recursion relation for the eigenvalues of Rg- 
If we call r\Q the only non-zero eigenvalue of Mq, we may write 

??g = XgAg = Xg[A^-i + ('7-1)Xg-i]- (68) 

From Eq. . one sees that Zq = r\Q^ since the trace is just this non-zero eigenvalue. 
Now, using Eq. it^ . we have 

XG = gXG-i> (69) 
from which we are finally led to the map 

^G = 9 {^l-x + (9 - l)x^-i] . (70) 
6. Thermodynamic functions 

If we take the Boltzmann constant /cb = 1, the free energy per monomer of the system 
may be written as 

T T 

fo = -^^nZc = -^Inijc- (71) 
Defining an auxiliary map, Kq ~ Xg/^Gj we have 

^ln(l + (g-l)if^_i, 2^ 

where 



/g = /g-1 -^ln{l + {q~ l)Kl_,) - :J In 2, (72) 



The recursive iteration of Eqs. H72|l and (|73|l . with the initial conditions fo — —e 
and Kq = exp(— /3e), leads to the free energy per monomer for any generation G, which 
converges to a well-defined free-energy in the thermodynamic limit, G — > cxd. Maps 
for additional thermodynamic functions, as the entropy and the specific heat, can be 
obtained by taking the derivative of Eqs. H72|l and (|73(l with respect to temperature. 
For the entropy per monomer, for example, we obtain 

' + 2^ (1 + (, - i)A-6-o + + i >..^- (») 

where 

dKo 2qKG-i dKa-i 



dT {l + iq-l)Kl_,f dT 



(75) 



As an example, in Figure 5 we show the results for the entropy and specific heat 
of a uniform model on the lattice with p = 2 and q = A. Numerical analysis shows 
there is a genuine singularity, associated with a divergence of the specific heat at a 
critical temperature. Note the interesting behavior of the system above the transition 
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Figure 5. Entropy and specific heat for a uniform model, in a lattice with p = 2 
and g = 4, calculated using the transfer- matrix technique. In this figure e/fcs = 1. 



temperature, with a constant entropy per monomer, and, consequently, zero specific 
heat. On physical grounds, this result should have been anticipated, since the polymers 
are completely unbound on the high-temperature phase, and the maximum amount of 
disorder is attained independently of temperature. The RG approach, however, does 
not yield such a global picture of the thermodynamics, which is possible only in the 
TM framework. 

To check the reliability of the method, and its compatibility with the RG results, 
we may compare the critical temperature and exponent it yields with those predicted 
by the renormalization-group calculation. For any value of q, the fixed point y* = q—1, 
with y — (exp/3e), gives the critical temperature 

with a critical exponent given by Eqs. Q and For q — i and e/ks = 1, we have 
Tc = 0.910239 ... and a = 0.290488 . . ..The numerical analysis of data in Figure 4 
leads to Tc — 0.910239 and a ~ 0.2905, which confirms the accuracy of the method. 

Now we remark that the models we are interested in include more than one 
interaction energy, depending on the position i along the path between the root 
sites. The basic steps of the TM scheme can be adapted in order to obtain the 
corresponding maps, although much attention has to be paid to all of the details. For 
instance, if we consider an aperiodic model with the presence of two distinct interaction 
energies, which are placed along the lattice according to the period-doubling sequence 
a — > a6, 6 ^ aa, the method requires the definition of two matrices and M^, the 
eigenvalues of which are 77^' and t/q^ . The maps for 77^' and 77^^ are written as 



(a) 



vl,^liV^G^-^ + iq-l)x%-i , (77) 



Aperiodic model for interacting polymers 



17 



and 



(78) 



Note that each one of these eigenvalues gives rise to a different partition function, 
associated with the choice of a or 5 as the initial letter to be iterated according to 
the inflation rule. We are always interested on sequences generated by the recursive 
application of the rule to the initial letter a. 

If the aperiodicity is induced by the four-letter Rudin-Shapiro sequence, a —^ ac, 
b dc, c — > ab, d — > db, the set of four maps for the eigenvalues are given by 

(79) 
(80) 
(81) 



and 



(a) 

Vg 


= q 


(a) (c) 

Vg-iVgU - 


Hq- 




(b) 

Vg 


= q 


(c) (d) 

Vg-iVg-i H 




^)x'g-i 


Vg 


= q 


(a) (6) 




^)x%-^ 


(d) 
Vg 


= q 


(6) (d) 

Vg-iVg-i - 


Hq- 


i)xS-i" 



(82) 



The free energy per monomer, along with its temperature derivatives, can be similarly 
defined, so that the singularity at the phase transition can be analysed directly. 



7. Discussions and results 

As discussed in Sec. 2, the RG analysis of the uniform model indicates the presence 
of a second-order phase transition for q > 2, with the specific-heat critical exponent 
given by Eq. @. In case of irrelevant aperiodicity, that is, when the diagonal fixed 
point has just one relevant eigenvalue, a is given by the same expression. For the 
period-3 sequence, however, the diagonal fixed point is completely unstable, and the 
two-cycle should be responsible for the critical behavior, as in the case of the spin 
models |2- In this case, the scaling analysis must be somewhat adapted to take into 
account that two renormalization-group iterations are needed in order that the system 
goes back to the vicinity of one of the two points that are part of the two-cycle ^U] . 
The result is simply that the specific-heat critical exponent is now given by 

Inp 



2-2 



InA' 



(83) 



where A is the leading eigenvalue of the linearized second-iterate of the RG recursion 
relations about any one of the two points of the attractor. Results of this analysis 
have already been given elsewhere and will not be repeated here. 

For the model with Rudin-Shapiro aperiodic interactions, we remarked above that 
there exist two non-diagonal fixed points together with the curve composed of two- 
cycles, and that the linearization of the second iterates of the recursion relations 
about any point on the two-cycles curve gives the same eigenvalues. For each 
q > 2 + \/2, we may therefore determine numerically which value of a Eq. (|83|l 
predicts, and then compare it with the direct analysis of the singularity which comes 
from the TM method. It is also possible to obtain a in the usual way, linearizing 
the recursion relations (first iterate) about the non-diagonal fixed points, using the 
leading eigenvalue that comes from the solution of Eq. H22() . The coincidence of the 
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Figure 6. Specific heat of the model with Rudin-Shapiro apcriodicity, in a lattice 
with p = 2 and 5 = 4, calculated using the transfer-matrix technique. In this figure 
is 100 times greater than the other energies. 

values is already an important hint of the correctness of scaling predictions, and we 
have indeed verified it for several choices of q. 

In Figure 6 we show the TM results for the specific heat in a lattice with q = 4, 
with a certain choice of interaction energies. The first interesting feature is the 
appearance of log-periodic oscillations in the low-temperature phase, as in spin models 

This is a natural consequence of the discrete scale- invariance of the aperiodic 
sequence (due to its self-similar character), which implies a natural rescaling factor 
in the renormalization group |12| . Several different values of the interaction energies 
must be separately analysed, and the net result is an exponent a = 0.252 ±0.08. Two 
points must be carefully stressed: first, that some choices for the energies (ca = Cd 
and £{, = Ec, for instance) give rise to an effectively periodic model, because of the 
symmetries of the Rudin-Shapiro rule, and should therefore be kept out of the analysis; 
second, a does not show any important dependence on the values (provided they are 
not of the form that makes the model effectively periodic, of course), which points to 
a true "aperiodic universality class" associated with the Rudin-Shapiro geometrical 
perturbation of the model.Now, for q — 4, the scaling result is a = 0.253692 . . ., in 
striking agreement with the TM value. The same scenario is present for several other 
values of q we have tested, what leads us to believe in the correctness of the methods. 

In conclusion, we have presented detailed renormalization-group and transfer- 
matrix calculations for a class of interacting polymer models on diamond-like 
hierarchical lattices, with aperiodically distributed coupling constants. Although 
straightforward, the exact renormalization-group analysis has revealed a surprising 
family of attractors in case of Rudin-Shapiro apcriodicity, and this prompted us to 
resort to the transfer-matrix formalism to check the scaling results. However, we had 
to develop a complete reformulation of this method, in order to be able to apply it to 
the polymer problem. The transfer-matrix calculations have confirmed the results of 
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the simple scaling analysis, but have also revealed peculiarities of the transition that 
were not accessible to the rcnormahzation-group study. What is most important is to 
notice that aperiodic perturbations may lead to new universality classes, adding up 
to the usual criteria of dimensionality and symmetry. In a sense, the breakdown of 
translation invariance may be relevant to the determination of new universal behaviors, 
and the particular way in which this invariance is broken must be taken into account. 
The introduction of disorder, for example, is a way of breaking translation invariance, 
and there are several instances in which its effects on critical behavior are well-known. 
Aperiodic distributions of couplings are just another way of accomplishing this, and, 
although more difficult to implement physically, they are amenable to more controlled 
calculations, such as those presented in this paper. 
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